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Abstract 

We investigate the quantum phase diagram of the exactly solved 
mixed spin-(^, 1) ladder via the thermodynamic Bethe ansatz (TBA). 
In the absence of a magnetic field the model exhibits three quan- 
tum phases associated with su(2), su(4) and su(6) symmetries. In 
the presence of a strong magnetic field, there is a third and full sat- 
uration magnetization plateaux within the strong antiferromagnetic 
rung coupling regime. Gapless and gapped phases appear in turn as 
the magnetic field increases. For weak rung coupling, the fractional 
magnetization plateau vanishs and exhibits new quantum phase tran- 
sitions. However, in the ferromagnetic coupling regime, the system 
does not have a third saturation magnetization plateau. The critical 
behaviour in the vicinity of the critical points is also derived system- 
atically using the TBA. 
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1 Introduction 



The field of exactly solved models in statistical mechanics has many signif- 
icant highlights. These include Elliott Lieb's pioneering work on the six- 
vertex model and his calculation of the residual entropy of square ice pQ. 
Over the ensuing years the six-vertex model and the related Heisenberg spin 
chain have been generalised in all manner of directions. Most recently at- 
tention has turned to the physics of quantum spin ladders, for which a num- 
ber of exactly solved models have been proposed. The underlying model is 
the spin- 1 Heisenberg ladder, which has been studied extensively This 
model consists of two Heisenberg chains coupled together with Heisenberg 
rung interactions forming a ladder-like structure. A number of ladder com- 
pounds have been synthesized, such as SrCu 2 3 [3], Cu2(C 5 H 12 N 2 )2Cl4 0j, 
(C 5 H 12 N) 2 CuBr 4 0, (5IAP) 2 CuBr 4 -2H 2 [6 , KCuCl 3 and TlCuCl 3 0. The 
experimental results reveal an interesting mix of low-temperature physics, in- 
cluding spin excitation gaps and magnetization plateaux. 

The theoretical investigation of the ladder compounds has also been cen- 
tred on a number of variants of the standard Heisenberg ladder j^j, includ- 
ing the addition of multi-body interactions [""], alternation and frustration 
[T""". IT""] . The ladder models have been studied by a variety of methods, 
for example, numerical [""J H2] , perturbation theory [T3J E] and the quantum 
transfer matrix algorithm [JJ"] . Unfortunately the Heisenberg ladder is not ex- 
actly solvable in the sense of the six- vertex model or related su(2) Heisenberg 
chain. However, a number of variants have been solved exactly by means of 
the Bethe ansatz; see, for example, Refs [TBI El CHI El EI11 EJ] • For arguably 
the simplest model, based on su(4), the critical behaviour derived from the 
thermodynamic Bethe ansatz (TBA) is seen to be consistent with the exist- 
ing experimental, numerical and perturbative results for the strong coupling 
ladder compounds [22] . This includes the spin excitation gap and the critical 
fields H c i and i7 c2 , which are in excellent agreement with the experimental 
values for the known strong coupling ladder compounds (5IAP) 2 CuBr4-2H 2 0, 
Cu 2 (C 5 H 12 N 2 ) 2 Cl 4 and (C 5 H 12 N) 2 CuBr 4 . 

On the other hand, the special interest in fractional magnetization plateaux 
[2*3] has inspired work on mixed spin chains [2"""] |2**5] EE1 EZ] > mixed spin lad- 
ders |2E] and various experimental compounds [2H1 HU31 EH]. In particular, 
the magnetic behaviour of a mixed spin-(|, 1) Heisenberg ladder has been in- 
vestigated by means of the density-matrix renormalization group technique 
[2*5] . It was concluded that for certain strong rung coupling magnetization 
plateaux exist at M z = 0.5 and at M z = 1, but with no plateaux for nega- 
tive (ferromagnetic) rung coupling. In fact the mixed spin ladder exhibits a 
richer phase diagram than the spin-| ladder. It appears that a comprehensive 
study of the mixed spin ladder, e.g., the prediction of the critical fields for 
different rung coupling and examination of the critical behavior, has not been 
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undertaken. It also remains to investigate the effect of non-equal Lande g- 
factors on the ladder legs, particularly given that g-factor anisotropy appears 
to affect the critical fields in the spin-orbital model [321 EH] • 

In this paper, we investigate the quantum phase diagram of the exactly 
solved mixed spin-(~,l) ladder with different ^-factors, via the TBA. 
We find that in the absence of a magnetic field the model exhibits three 
quantum phases associated with su(2), sw(4) and su(6) symmetries. In the 
presence of a strong magnetic field h, two magnetization plateaux appear in 
the strong antiferromagnetic rung coupling regime. The fractional magneti- 
zation plateau M z = \g s ^B+\{.gs—9t)^B corresponding to the fully-polarized 
doublet rung state, opens at the critical field H c i and vanishes at the crit- 
ical field H c2 . The second plateau M z = (~g t + g s )fiB, corresponding to a 
fully-polarized quadruplet rung state, opens only at a very strong magnetic 
field H c3 . For weak antiferromagnetic rung coupling, the fractional plateau is 
closed such that three different kinds of quantum phase transition occur. For 
ferromagnetic rung coupling only the full saturation magnetization plateau 
exists. The critical behavior for the different quantum phase transitions is 
systematically derived by using the TBA. 

The paper is organized as follows. In section [21 we discuss the exactly 
solved mixed spin-(|, 1) ladder model with different Lande factors. The exact 
solution is given via an appropriate choice of rung basis. Section El is devoted 
to the investigation of the quantum phase diagram. In Section 0] we investi- 
gate the magnetization plateaux in the presence of a strong magnetic field. 
The critical fields characterizing the different quantum phase transitions are 
given explicitly. A summary of our main results and conclusions are given in 
Section El 



2 The exactly solved mixed spin- (|, 1) ladder 
model 

The Hamiltonian of the exactly solved spin-(~, 1) ladder model, based on the 
su(6) symmetry, reads [TH] 



L L L 

H = J {] H leg + J±J2 ■ 4' - 9t^h J2 T J~ ^ S] } (1) 

j=i j=i j=i 



7=1 ^ ' 



(2) 



— » — » -, 

Here Tj and Sj are the standard spin-^ and spin-1 operators acting on site 
j of the upper and lower legs, respectively (see Fig.l), J\\ and J± are the 
intrachain and interchain coupling constants, h is the magnetic field, //b 
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Figure 1: The mixed spin ladder: Jn and J± are the intrachain and interchain 
couplings; g t and g s are the Lande factors along each leg. 

is the Bohr magneton and gt and g s are the Lande factors along each leg. 
Periodic boundary conditions are imposed with L the number of rungs. 

The rung term in (Q) breaks the sit (6) symmetry of ifi eg into su(4) ©w(2) 
symmetry. This symmetry is in turn broken under the magnetic field. The 
physical properties and the critical behaviour of the model are determined 
by the competition between the rung and leg coupling constants and the 
magnetic field h. Essentially, H\ eg is the permutation operator corresponding 
to the sit (6) algebra symmetry. If we change the canonical basis e^ 13 <g> ej s 
of V\ (g> V 2 into rung quadruplet and doublet states (Clebsch-Gordon decom- 
position), the six- dimensional space splits into the direct sum of quadruplets 
and doublets with the basis 

n> = ^(il-D-^J)). l2> = ^ (l -M>- ^lo, -t>), 

|3> = |4) = ^(|l,-|) + v / 2|0,I)), 

|5> = 75(1-1,!) + V^iO, -§)) , |6> = I - 1, -i>, (3) 

where the states |1), \2) form the doublet and the remaining states form the 
quadruplet. The projectors onto the doublet and quadruplet subspace are 
given by 

P d = _f(f •£-§), P q = |(f -S + l). (4) 

It follows that the rung interaction term can be accommodated into an su(Q) 
invariant Heisenberg chain by embedding the doublet rung states through an 
appropriate chemical potential term. The leg and rung part of the Hamilto- 
nian ([TJ can be derived from the relation 

H = Jii-f- In t(v)\ v=0 + E ±+h + const (5) 
dv 

associated with the quantum transfer matrix r(v) = troT(u). The energy 
E± + h arising from the rung interaction and magnetic field terms is given 
further below. Here T(u) denotes the monodromy matrix given by 

T(v) = R , L (v)R 0jL ^(v) . ..Ro, 2 (v)Ro,i(v) (6) 
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associated with the su(6) quantum i?-matrix. 

Now consider the effect of the magnetic field. Although the magnetic 
field preserves the integrability of the leg part of the Hamiltonian, the dif- 
ferent ^-factors on each leg break the doublet /quadruplet basis J3J) for the 
Hamiltonian (^Q). Fortunately, we can still find another basis, 

V>§ = |1,|> S ^-§ = 1- 1,-1), 

to diagonalize the rung and magnetization terms simultaneously. The quan- 
tities are given by 



y&> = aV2[(g s - g t )h' + a] ± y/l + i (</./*' - </ t /i' + a) 2 , (8) 

where a = ±| and ft, = //b/i/Jl. We notice that if g s = g> t , the basis states 
^ and reduce to the doublet, with the other states reducing to the 

2 2 

quadruplet. With regard to the total spin of the multiplets we can still call 

the states (ip3 , , ip_s) quadruplets and (ip^~~\ doublets. 

2 2 2 2 2 2 

It is well established that the Hamiltonian (JTJ) can be diagonalized via the 
algebraic Bethe ansatz. In this procedure it is important to note that the 
leg content of the Hamiltonian, H\ eg , is not altered under the change of basis 
order between the quadruplet and the doublet states, however the rung and 
magnetic field terms are altered by these changes. We note also that for the 
ladder Hamiltonian (JTJ) the doublet rung state is energetically favoured for 
J± > 0, whereas the quadruplet rung state is favoured for J± < 0. This is 
the reason for choosing the doublet state as reference state for J± > 0, while 
a quadruplet state is choosen as reference state for J± < 0. As the magnetic 
field is turned on, the energy levels of each multiplet component split. The 
basis order is therefore chosen in accordance with their energy levels. 

The resulting Bethe ansatz equations are well known jSlj and consist of 
a set of five coupled equations depending on five flavours, v^ k \ k = 1, . . . , 5. 
The Bethe ansatz equations 

n v j v i + 2 = TT v j v l + 1 TT v j v l 2 / Q N 
(k) (k-l) i 11 (fc) (fc) .11 (fc) (fc+1) . j V 3 ) 
i=l V ) V i 2 1=1 V ) ~ V l - 1 1=1 V ) ~ V l +1 

can be derived from the nested algebraic Bethe ansatz. In the above, k = 
1, . . . , 5 and j = 1, . . . , M k and the conventions v ^ = vf^ = 0, M 6 = apply. 
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In real ladder compounds the difference between the ^-factors along each leg 
is small (this is not always true for spin-orbital models |22])- Thus from now 
on we treat the term g s — g t as a small quantity. After some algebra, the 
eigenspectrum is obtained from relation (J3J) as 



E 



Mi 



i=l V t + 4 



+ E 



±+h, 



(10) 



where the energy contribution from the rung interaction and the magnetic 
field terms is given by 



E 



±+h 



\Jx- \g s ^Bh - ^= J±\/l + l(g s h' - g t ti + 



1)2 
2 - 1 



(-) 

2 



+ 



+ [|J± - (5</t + <7«)^] 
-\J±. ~ \g s ^sh + -^J±\j 1 + I - 9th' + 



1)2 

2^ 



1)2 

2/ 



(+) 



+ 



+ (f#t + # s )/wi]N-§ 



1)2 
2^ 



iV 



(+) 



(-) 



-(l^ + g s )^BhNs - [\g t n B h - l(g s - gt)^Bh]N[ 
+[\g s V>Bh - \{g 3 - gt)i^Bh}N { V 
+{\gt + g s )HBhN_3 + const. 



(+) 
i 

2 



'ID 



Here the iV's are the numbers of the corresponding states. In the thermody- 
namic limit, the Bethe ansatz equations © allow the string solution 

- {k)n - " ^; i + |i(n + l-2i), 



(12) 



where j = 1, n, a a = 1, N n and , k — 1, 5, are the positions of 

the center of the strings of flavour k. The number of n-strings, N n a \ satisfies 
the relation = *Yl, n nN n k \ On taking the thermodynamic limit, the 

Bethe ansatz equations become 



rn 



rn 



P {2) 
rm 5 



a n ^ A nm * pU + ^ 0>nm * 
m m 



(13) 
(14) 



where k = 2, . . . , 5, and the symbol * denotes convolution. Pn\v) and 
pn^ h {v) with k = 1, ... ,5 are the densities of roots and holes for the five 
flavours. We have adopted the standard notations 

A nm {\) = S(X)S nm ~ V u nmJ u '\n—m\ (A) + (A) (15) 

Min(n,m)-1 

+2 ^2 a\n-m\+2l(X), 

1=1 
Min(ra,m) 

(A), (16) 

z=i 

with a n (X) = ^^t^. 

In order to find the equilibrium state of the system at a fixed temperature 
T and external magnetic field h > 0, we minimize the free energy F = E — 
TS — hM z with respect to the densities and then obtain the TBA equations 
in the form 

e (k) oo £ (fc) 

e [ k) = ^f ) + Ta 2 *ln(l + e-^)+T(ao + a 2 )^a m *ln(l + e-^ i ) 

m=l 

°° / (fc-i) (fc+l) \ 

-T ^ a m * ( ln(l + e~^) + ln(l + e~^) J , (17) 

m=l ^ ' 

ei fc) = ^ ) +Ta 1 *ln(l + e !2 ^ i +Ta 2 *ln(l + e-^)) 

+T(a + a 2 ) ^a m _ n *ln(l + e-^) (18) 

m>n 

°° / e (*-l) e (fc+l) \ 

-T^a m _ n+1 * hn(l + e-^) + ln(l + e-^)J , n > 2. 

m>n ^ ' 

Here pi fc)/l (A)/pi fc) (A) := exp(ei fc) (A)/T) with fc = 1, ... ,5 and ei 0) (A) = 
£n\X) = is assumed. The dressed energies play the role of excita- 
tion energies measured from the Fermi level for each flavour. The driving 
terms in the antiferromagnetic rung coupling regime for a weak magnetic 
field, h < -, — p3 — — ™ — , are given by 

= ~ J \\ V 2 l + 1 + \9s + \{9s ~ 9t)]pBh, 

gf ] = P± - [\gt + \g, + l(g s - #t)W^, 

g ( i } = [k(gt + gs) + l(g s -gt)}PBh, (19) 

5-i 4) = [& - §(& - 9t)]PBh, 

g? = [\{gt + 9s) + l(g s - g t )}pBh. 
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The higher driving terms are given by gk = n[g s + \{g a ~ fiO]/AB^ an d 
grf* = ng[ k ^ for k > 1. Consequently, the free energy is given by 



= -\g s ^ B h -\{g s - gt)^ B h -T a «( A ) ln d + e~^)dA.(20) 



It is worth mentioning that the driving terms vary for different choices of 
the basis order. The TBA equations (|T7j) and (f*H?j) provide a clear physical 
picture of the groundstate and make the thermodynamic properties, such as 
the free energy, magnetization and susceptibility accessible. 

3 Quantum phase diagram 

In the low temperature limit, the states with positive dressed energy are 
empty. The rapidities with negative dressed energy correspond to occupied 
states. The zeros of the dressed energies define the Fermi energies. As usual, 
we decompose into positive and negative parts, ej? = + e« ~ , with 
only the negative dressed energies contributing to the groundstate energy. 
Analysis of equations (fTTj) and (fTU|) in the limit T — > reveals that the roots 
are all real for the groundstate, corresponding to n — 1. All dressed energies 
e„. with n > 2 correspond to excitations. Under this circumstance, we 
see that all energy bands are completely filled in the absence of an external 
magnetic field and rung interactions. 

In order to derive the groundstate properties, we first consider the anti- 
ferromagnetic regime J± > 0, where the doublet component ip[ ■* is chosen 



as the reference state with a basis order ,ip_^,ip3,ip[ \ip_l 1 ip_-±). The 

2 2 2 2 2 2 

groundstate TBA equations then read 



It is clear that without a magnetic field the su(6) multiplet levels split due 
to the rung coupling. If J± is very large (the limit of strong rung coupling), 
the whole quadruplet state (ip3 , ^i^i' ^-2 ) i s gapfull, i.e. > for 

2 2 2 2 

k > 2. This means that the quadruplet is not involved in the groundstate - 
the groundstate consists of doublet states with massless excitation. Solving 
the TBA equation (|2~T|) . we find that the quadruplet excitation gap is given 
by Ai = | J± — 2J|| In 2. Thus if J± becomes larger than the critical value 
J+ = |J|| In 2, there is a quantum phase transition from the Luttinger liquid 
su(A) © su(2) phase into the su(2) phase. However, in the presence of a 
magnetic field this critical point is not stable. In this case the dressed energy 





M- +ai * e (2)- 

•(*)- + 0l * p-D- + e <*+l)-] , 



(21) 
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Figure 2: The doublet polarized state forms a rung trimerized groundstate. 

levels become completely split by the magnetic field h. If the rung coupling 
J_i_ is large enough so that the driving term in remains positive, 
the quadruplet state could be gapfull and the groundstate would still be the 
doublet. However, as the magnetic field increases, the doublet component 
ip^i gradually shifts out of the groundstate as the Fermi surface of the dressed 

2 

energy e^ 1 ) lifts. Subsequently, if the magnetic field h is larger than the 
critical point H&, the reference state becomes a true physical state such that 
the strong rung coupling forms a rung trimerized groundstate (see Fig. 2). In 
other words, the doublet component ip^ forms a ferromagnetic groundstate. 

2 

The critical field H c \ is given by 

H cl = AJ ll /[g s + l(g s ~g t )]fx B . (22) 

The critical point H c i indicates a quantum phase transition from a gap- 
less magnetic phase into a ferromagnetic phase with gap A 2 = (\g s + \{g s — 
g t )]fiB{h — H c i). It is worth noting that in this ferromagnetic phase a magne- 
tization plateau M z = \g s ^B + \{g s — 9t)^B opens. The necessary condition 
for this plateau to exist is 

r ^ j+f 8J n [\gt + lgs + \(g s -gt)} 

J ± > J c = F ~T7 yj • { I6 > 

3 [9s + sKQs- gt)\ 

In the critical phase h < H c i, the TBA equations with a very large or a very 
small Fermi boundary can be solved analytically. If the magnetic field is very 
small, i.e. /id, the Fermi boundary of the dressed energy is very large. 
The energy potential satisfies the Wiener-Hopf type equation 

IT f~ B,oo 

e(1) ( A ) = - J \\—r^ + \\9> + l(9s - 9t)]n B h + / G(X — k)* e^(k)dk. 

COStlTlA J-oo,B 

(24) 

Here the function G(X) is defined via 

G(X) = — / —— -e-'^du. 

v 1 2tt 2cosha;/2 
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Using the standard Wiener-Hopf technique, we find that the Fermi boundary 
satisfies the relation e~ BlT = [g s + |(g s — g t )]fi B ha-(0) /4J\\ira + (m), where the 
decomposition functions are given by 

M<-)=M-) = ^(^)"*/T(i-|). P5) 

The Fermi boundary decreases monotonically with increasing magnetic field. 
Correspondingly, the free energy is given by 

» -J„[*(l) - *(!)] - JLfo. + l(g s - g t )] 2 n 2 B h\ (26) 

which suggests a susceptibility of % 4^2 [fl^ + |(<7s — <7*)]Vbj indicating 
an su(2) critical phase. In the above T(z) and \&(a) are the gamma and 
diagamma functions, respectively. 

On the other hand, if the magnetic field h tends to the critical point H c i, 
the Fermi boundary Q of the dressed energy e^ 1 -* is very small, say Q <C 1 for 
H cl — h<^l. Under this circumstance, the free energy is given by 

F(y °' — ~ -y s VB - l(g s ~ gtjUBh - — [g s + l(g s - g t )}/j, B (H cl - h), (27) 

Li 7T 

where Q a/ (H c1 — h)/AH cl . Thus the susceptibility 

x ^^^M iHcl -h)-i (28) 

indicates the singular behavior of the transition from the gapless phase into 
the gapped phase. In addition, from the Bethe ansatz equations and the 
relation 



M z 



rQ 

\g s HB + \{g s - gt)nB-[g s + \{g s - gt)\nB I p?(\)d\, (29) 

J-Q 

the magnetization per site M z in the vicinity of H cl follows as 

40 20 
M z = \g sl i B + \{g s - g t )n B - -% + \{g s - g t )]fi B (l - -^). (30) 

71 71 

Apparently, as h—>H c i the magnetization M z tends to the plateau value 
\g s ^B + |(<?s — gt)^B- We shall discuss the magnetization plateaux as well 
as the quantum phase transitions in the next section. 

For the ferromagnetic regime J± < 0, the rung quadruplet component ibz 

2 

is chosen as the reference state with (ip3 , V^ 3 - > Vd as the 

2 2 2 2 2 2 
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Figure 3: The quadruplet fully-polarized state forms a trimerized ferromag- 
netic groundstate. 



order of the basis. Thus the driving terms are given by 



9i 



(i) 



-J\ 



1 



Ugt + 9s) + Ug a - 9t)]f*Bh, 



9i 



(2) 



9i 



(3) 



9i 



(4) 



9i 



(5) 



v 2 + \ 
[9s - \{g s -gt)]^ B K 
[\{gt + g s ) + \{g s - 9t)]^Bh, 
-\J± - [\gt + \g s + \{g s - g t )]^Bh, 
[g s + l(g s - gt)\v>Bh. 



(31) 



In the absence of a magnetic field, the quadruplet and doublet states are 
degenerate. If the rung coupling becomes negative enough the doublet state 
is completely gapfull, and the Fermi boundaries of the quadruplet states are 
infinity. Using Fourier transforms, we find that the doublet does not exist in 
the groundstate for J± < J~ = — |J||. Again this critical point is not stable 
if the magnetic field is applied. We can see that the Fermi surfaces of the 
dressed energies e^ k \ k = 1,2,3 lift, while the Fermi surfaces of the doublet 
sink. If the magnetic field is strong enough, i.e. if 



h>h c = 4J||/[§(# t + g s ) + \(g s - g t )}fi B , 



(32) 



the reference state ip3 becomes a true physical state. Thus the groundstate 
is a fully-polarized (trimer-like) ferromagnetic state (see Fig. 3). Note that 
the rung coupling J± must be less than a critical value JJ F given by 



8Jn 



[\gt 



+ \g s + l(g s -gt)\ 



(33) 



4 Magnetization plateaux 

Magnetization plateaux are one of the most interesting phenomena in the lad- 
der compounds. For example, fractional magnetization plateaux have been 
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found in Shastry-Sutherland systems j^Zj- Theoretical studies and numeri- 
cal results suggest that magnetization plateaux exist in the mixed spin-(~, 1) 
chains |2*H 1251 12*0] and the mixed ladder [2E]- From the analysis of the critical 
points of the solvable model in the last section, we found that gapped or gap- 
less states appear in turn as the external magnetic field increases. For very 
strong rung coupling, i.e. J± ^> J^ F , the two-component massless quantum 
magnetic phase lies in the regime h<H c i. The ferromagnetic phase appears 
for a magnetic field h > H c i and the component ip[ ^ becomes a physical 

2 

ferromagnetic groundstate. On the other hand, the magnetic field can bring 
the state i[)3 close to the groundstate. Eventually it becomes involved in the 
groundstate when the magnetic field is strong enough. From the expression 
(fTTj) . we see that if h > 3J±/ (g t + 3g s + ~(g s — g t ))jiB, the state ipz becomes 

an energetically lower lying state than In this case it is convenient to 

2 

reorder the basis as ipi , > ip-i )• Here the TBA driving 

2 2 2 2 — 2 2 

terms are given by 

9i ] = ~J\\ i +\J±~\\{9t + 9s) -\{9s- 9t)}HBh, 



V 



4 



9i = + [\9t + %9a + \(jg* - 9t)]^Bh, 



9? = P±-gsVBh, (34) 



9 { i ] = [9s - Ug s - 9t)]HBh 



9i ] = [\{9t + 9s) + l(g a - 9t)}^Bh. 

Analysing the TBA with these driving terms, we find that the ferromag- 
netic groundstate can be maintained only in the regime H c i < h < H c2 , 
where the critical point H C 2 is given by 

3 J \ — 8Jii 
[{9t + 9s) - ^{g s -gt)\^B 

Beyond the critical field H c2 , the plateau M z = \g s HB + \{.g s —gt)^B vanishes. 
In the vicinity of H C 2, i.e. h — H C 2 <C 1, the Fermi boundary is very small. 
After a similar calculation the Fermi point is found to be 



^ ,[\(9t + 9s) -\{9s- 9t)]^B{h- H c2 ) 

Q~\l igT\\ • (36) 



The susceptibility is given by 



_ 3 [|(& + 9t)^B ~\{9s- 9t)liB] 2 (37) 
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which indicates the singular behavior of the transition of the gapped phase 
into the gapless phase. The magnetization shows the square root field de- 
pendent behaviour 

40 20 

M z = [\g s + \(9s- gt)\»B + —[\{gt + g s ) -\{g.- <?*)W(l - — )• (38) 

7T 7T 

If the magnetic field is further increased to h > 3J±/[(g t + g s ) — |(# s — 
g t )]fj,B, the state if) a becomes the lowest lying state. Thus in this regime, 

it is reasonable to choose the basis order as (ip3,ip{ \ ^ + , > V'-i)- 

2 2 2 2 2 2 

Subsequently, the driving terms are given by 

g[ 1] = -J\\ 2 1 i - f J± + [\{gt + g a ) - \{g s - gt)W/t, 

g[ 2) = ~J± + ~(g s -g t )]vBh, 

g? = -lJ± + g s ^ B h, (39) 

g{ = \Jx-\{g s - gt)t*Bh, 

gf = [§(& + &) + l(9s -9t)]PBh. 

As the magnetic field increases beyond the critical point H C 2, the groundstate 
becomes a mixture of the doublet and quadruplet states. Strictly speaking, 
the doublet component ib^ and the quadruplet component ips compete for 

2 2 

the groundstate. Other components of the multiplets are gapfull by virtue 
of both the rung coupling and magnetic field. As the magnetic field becomes 
stronger, the probability of the quadruplet component ip3 becomes higher. 

The inflection point at hjp = 3J±/[(g t + g s ) — \ {g s — gt)]^B indicates an equal 

probability between the components ip3 and which can be seen clearly 

2 2 

from the magnetization curve in Fig. 4. Using the TBA equations with the 
driving term (|39jl . we find that for a sufficiently large magnetic field h > H c3 
the groundstate becomes fully-polarized with a full saturation magnetization 
plateau at M z = (^g t + g s )^B- The critical point is given by 

3 J ± + 8 J\\ , s 

[{gt + g s ) - ^{g s - gt)\fi B 

Analogously, we find the singular behavior in the vicinity of the the critical 
point H c3 . The susceptibility is given by 



3[|(& + gt)»B - \{g s -gt)^ B ] 2 

Any/JwlHa - h) 

which indicates the nature of the singular behavior of the transition between 
the gapless and gapped phases. The magnetization also exhibits the square 
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Figure 4: The magnetization versus magnetic field h in the strong antiferro- 
magnetic rung coupling regime. The magnetic moment is normalized from 
the magnetization via M z /^(g s + gt)/J-B- The coupling constants used are 
Ji = 6.0K, J|| = OAK and we take (/-factor values g s = 2.22, g t = 2.09 
with hb = 0.672K/T. From the TBA we predict that the one third satu- 
ration magnetization plateau opens only if J± > J+ F as given in (J23|) . The 
indicated critical fields H cl « 1.052 T, H c2 « 5.162 T and H c3 « 7.394 T 
predicted by the TBA coincide with the numerically estimated values. The 
inflection point at h = hyp ~ 6.278 T and M ~ 1 indicates a point of equal 
probability for the states ip3 and ?/>i ^ ■ 

2 2 



root field dependent behaviour 

M z = [\g t + g^B - ^-[\{g t + g.) - \{g s - g t )]» B (l - % (42) 
where the Fermi point Q is 



71 Z " " 71 




\{gt + 9s) - \{g s - gt)] Mgg - /Q 



16J|| 

The magnetization increases almost linearly between the critical fields i7 C 2 
and if c3 . 

We have obtained the whole magnetization curve by numerically solving 
the TBA equations in the different phases (see FigJU). A third and full satu- 
ration magnetization plateaux is observed. On the other hand, as mentioned 
above, the first magnetization plateau at M z = \g s ^B + \{g s — 9t)^B de- 
pends mainly on the rung coupling. If J± < J+ F this plateau disappears. 
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1 r 




Magnetic Field H(T) 

Figure 5: The magnetization versus magnetic field h in the weak antifer- 
romagnetic rung coupling regime. The numerical values are the same as 
for the previous figure, but now using the smaller rung coupling constant 
J± = 1.3K < J+ F 2.07K. In this case the fractional magnetization 
plateau vanishes. The TBA critical fields coincide again with the numeri- 
cally estimated values. The inset shows an enlargement of the magnetization 
between H c \ and H c2 . 
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J+ F is a necessary condition for the existence of the one third saturation 
magnetization plateau. This conclusion is reasonable because the leg part 
of the Hamiltonian (0) is stronger than the rung part, due to the multi- 
spin interaction terms in eq.([2|). Therefore, the rung coupling must be large 
enough in order to dominate the low temperature groundstate. If the rung 
coupling fulfills J+ <J±< J+ F the field-induced fractional plateau vanishes. 
In Fig. Owe show the full numerical magnetization curve for weak rung cou- 
pling J± = 1.3K. We clearly see that the fractional magnetization plateau is 
closed. From the TBA analysis the first critical point H c \ lies in the interval 

3(Jjl - J+) ^ u ^ 3J ± 



[g t + 3g s + (g s - g t )/3}nB [gt + 3g a + {g s - 0t)/3]// B ' 

which for the given parameter set (see figure captions) evaluates to 0.47 < 
H c \ <0.66 T. Indeed, from the numerical results we find H c \ « 0.605 T. This 
implies that for H < 0.605T the groundstate is the doublet spin liquid phase. 
A magnetic field beyond this point allows the quadruplet component ipz to 
be involved in the groundstate. Therefore the critical point H c i indicates a 
quantum phase transition from a two-state phase into a three-state phase. 
Hence for h> 0.605 T, two Fermi seas, and e^ 2 \ lie in the groundstate. We 
can see that the probability of the component ijjz to be in the groundstate 
increases as the magnetic field increases. Meanwhile the doublet state ip^i 

~ 2 

is quickly driven out of the groundstate at the critical point H c2 , where 

the phase transition from the three-state into the two-state phase transition 

occurs. From Eq. 1)11)1. we find the middle point between if cl and H c2 is 

approximately at 0.66 T, which is consistent with the numerical curve in 

Fig. El Thus from the theory we predict the critical field to be H C 2 = 

6J±/[(9t + gs) — ^(gs — gt)]^B — H c \ ~ 0.715 T. This again is in agreement 

with the numerical value of 0.713 T. In the region H C 2 < h < H c ^ the two 

components ?/>7 and ipa compete to be in the groundstate. If the magnetic 

2 2 

field is strong enough, so that H > H c3 where H c3 is given by (J3UJ), the 

reference state ipa becomes a true physical groundstate. 



5 Conclusion 

We have investigated the phase diagram of the exactly solved mixed spin- 
(|, 1) ladder model in both the absence and presence of an external magnetic 
field using the TBA. It has been shown that in the strong antiferromagnetic 
rung coupling regime there exists a third and full saturation magnetization 
plateaux. A Luttinger liquid magnetic phase exists in the regime h < H c %, 
which corresponds to the doublet su(2) phase. The magnetic groundstate 
consisting of two components ipi and ip3 lies in the regime <h< H c3 . The 

2 2 
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ferromagnetic ground state with a third saturation magnetization plateau ap- 
pears in the regime H c \ <h< H C 2. The full saturation magnetization plateau 
opens at h > H c3 . The gapped or gapless states appear in turn as the mag- 
netic field is increased. The weak rung coupling regime exhibits three dif- 
ferent phase transitions, which involve two- and three-state quantum phase 
transitions. The fractional magnetization plateau vanishes. The model does 
not exihibit a third saturation magnetization plateau in the strong ferromag- 
netic rung coupling regime. We have also investigated the singular behaviour 
in the vicinity of the critical points via the solutions of the TBA equations. 
As the contributions from the leg interaction to the groundstate energy are 
very small in the strong rung coupling regime we believe that the solvable 
model ((TJ) is well suited to describe the physics of real mixed spin-(~, 1) lad- 
der compounds with a relatively large rung coupling constant. However, such 
mixed spin ladder compounds are yet to be found. 

One compound that we are aware of is the organic ferrimagnet PNNBNO 
p?T] , which has been recognized as a ladder compound with alternating spin- 
| and spin-1 units, i.e. as two coupled alternating mixed spin chains. The 
strong interchain (rung) coupling suggests that PNNBNO can be effectively 
identified as a mixed spin-(i, 1) ladder model via the Hamiltonian Q in the 
high temperature limit (T > 50 K). In this case the dominant rung interac- 
tion forms an effective one- dimensional spin-| antiferromagnetic chain. Of 
course for low temperatures the ferrimagnetic correlations between the spin-1 
and spin-| units should appear. Nevertheless we believe that the mixed spin- 
(~, 1) ladder model discussed here can at least describe the high temperature 
properties of PNNBNO HH]. 
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